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Abstract. Analyzing the archival data from the Rossi X-ray Timing Explorer (RXTE), we study the power density spectra 
(PDS) of Cygnus X-1 from 1996 to 2003 in the frequency range of 0.01 - 25 Hz. Using a model consisting of one or two 
Lorentzians and/or an exponentially cut-off power-law, we are able to achieve a good fit to the PDS during the observations. 
With our model we are also able to track the evolution of the Lorentzian components through all spectral states of the source. 
We confirm the relation between characteristic frequencies seen both in black hole candidate and neutron star sources, and 
show the changes in this relation during the transitional and soft states of the source. The connection between the Lorentzian 
components is investigated by analyzing similarities and differences in their behavior. We find that the spectral state of the source 
can be uniquely determined from the parameters of the these components. The parameter correlations can all be described by 
continuous functions, which differ between components. We discuss our results in the context of relativistic precession model 
for the accretion disk, and show a remarkable agreement between the model prediction and the data in the hard state. We 
estimate a value for the specific angular momentum of a, = 0.49 (-0.57) in the case of prograde (retrograde) rotation and an 
estimate for the inner radius of 22 to 50 (25 to 55) gravitational radii. Additional assumptions are required to explain the soft 
state data, and attempting to invoke rotational reversal for state transitions shows that it is insufficient to explain the differences 
between the hard and soft state data. 
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1. Introduction 

In many X-ray binary systems, the existence of distinct spec- 
tral states has been well documented. Often quoted as a pro- 
totype black hole system, Cygnus X-1 has been extensively 
observed and modeled. The source exhibits two main spectral 
states, usually referred to as hard and soft, with a brief inter- 
mediate state (sometimes interpr eted also as a very high state, 
see e.g. discussion in Zdziarski & Gierliriski 2004, and ref- 
erences therein) during the transitions. Several models have 
been proposed to explain the spectral states and transitions, 
with the two main components being a geometrically thin, op- 
tically thick accretion disk and a hot inner flow or corona. 
Seed photons from the disk are bel ieved to gain ener gy in 
the corona through Comptonization ( Shapi ro et alJ ll976). The 
models vary in the geometry and properties of mainly the 
corona/comptonizing region (for a review on coro nal mod- 
els, e specially in regard to timing characteristics, see lPoutanenl 
l200ll) . However, changes in the inner radius of the accretion 
disk are a common source of variability in many models (e.g . 
Poutanen et ail 1 19971 lEsin et aljfl998t IChurazov etaD 1200 it 
Zdziarski et alJ2002h . 
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The changes in spectral state are evident also in the power 
density spectrum (PDS). In the hard state the PDS is charac- 
terized by a flat-topped component up to ~ 0.2 Hz, where 
it breaks to a f -1 slope that steepens at a few Hz (e.g. 
iBelloni & HasingeJl99rilNowak et al1ll999l) . In the soft state 
the PDS is dominated b y a f ~ l component that steepens around 
10 Hz (e.g. lCui et al.ll997al) . The 1996 state transition also re- 
vealed an intermediate PDS, showing a ~ / slope at lower 
frequencies, a flatter component around 0.3 - 3 Hz, and above 
that similar to the hard state PDS. In many observations there 
is evidence for a quasi-periodic o scillation (QPO) at a few Hz 
(IBelloni et alll996UCui et alll997rllCui Il999h. 

Previous studies of the PDS of Cyg X-1 have shown 
correlations between temporal and spectral components (e.g. 
iGilfanov et alJll999i) and also between different temporal fea- 
tures. Correlatio ns between features in the hard s tate PDS were 
first reported by Wiin ands & van der Klisl dl999b . To resolve 
the origin of the tem poral features, it is vital to understand their 
long term behavior. IPottschmidt et alJ J20031 hereafter P03) 
recently presented a large scale study of the hard state PDS 
and also included data from 'failed state transitions'. However, 
no study has been presented where the PDS for all spectral 
states of Cyg X-1 are consistently modeled. To this end, we 
conducted a systematic analysis of all the available data from 
the Rossi X-Ray Timing Explorer (RXTE) satellite observations 
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of Cyg X-l in the HEASARC public archives. Using a model 
consisting of a power-law and two Lorentzian profiles, we were 
able to fit the majority of PDS, and to follow the components 
from hard state through the transitions and back. This allowed 
us to expand on previous results and monitor the evolution of 
the components. 

Our aim in this paper is to provide a consistent description 
of the PDS of Cyg X-l in all states, and through this study the 
evolution of the PDS throughout the time it has been observed 
by RXTE. We begin with a brief description of the available 
data, the computation of lightcurves and power spectra, and a 
presentation of our models in Sect. [5] This is followed by a 
presentation of our results, where we show the behavior of the 
different model components (Sect. [5}. In Sect.@]we follow with 
a discussion of the components, as well as putting our results 
in the context of physical models. Finally we summarize our 
findings in Sect. [5] 

2. Observations and Data Analysis 

2.1. The archival data 

Since its launch on December 30, 1995, and until the end of 
April 2004, the RXTE satellite has observed Cyg X-l in more 
than 750 pointed observations. The satellite has an orbital pe- 
riod of about 5.8 ks, but due to its low circular orbit, typical on- 
source time in a pointed observation is 3 - 4 ks. The satellite 
also carries an All-Sky Monitor (ASM) instrument on board, 
which continuously scans the sky, covering 80% every orbit. 

In this study we have used all public pointed obser- 
vation s, with data from the Proportional Counter Array 
(PCA. I.Tahoda et alJll996h . The PCA consists of five identical 
Proportional Counter Units (PCUs), and is sensitive in the 2-60 
ke V energy range. The observations span the time from January 
8, 1996 to April 3, 2004. The data mode used was the 'Generic 
Binned' mode. 1 For the count (hardness) ratios (see Sect. 12.41 
we used data in the Standard2 configuration. Due to the dif- 
ferent configurations and epochs for the observations, it was 
not always possible to obtain an exact match between energy 
ranges. We found that choosing the energy range of ~ 2-9 
keV allowed us to maintain the same range fairly well through- 
out the observations, with any deviations appearing in the upper 
boundary. We note that the shape of the PDS of Cyg X-l is not 
very sensitive to the wid th of the energy band in this region 
JNowaketal]ll999llP03l so the slight variations in the upper 
boundary should not affect our results. For the hardness ratios, 
the use of Standard2 data enabled us to keep the same energy 
ranges throughout the epochs. The channels and corresponding 
energy ranges are presented in Tabled 

We used the standard RXTE analysis software FTOOLS, 
Version 5.2. Lightcurves and X-ray spectra were extracted us- 
ing standard screening criteria: a source elevation > 10°, a 
pointing offset < 0!01 and a South Atlantic Anomaly exclu- 
sion time of 30 minutes. For the PDS, the lightcurves were 
extracted with a time resolution of 10 ms (for a few observa- 
tions a time resolution of 16 ms was the highest possible, due to 

1 For one observation run in epoch 1, PI 0240, only data from the 
single-bit mode was available, with the energy range ~ 2 — 9.5 keV. 



Table 1. Channel and energy ranges used in the count ratios 
and PDS. 





Count ratios 


PDS 




Lower band 


Upper band 


Upper channel" 


PCA Epoch 1 








channel 


0-10 


30-60 


35 


energy (keV) 


-2-4 


- 9-20 


9.5 


PCA Epoch 2 








channel 


0-6 


25-54 


29 


energy (keV) 


-2-4 


- 9-20 


9.3 


PCA Epoch 3 








channel 


0-6 


21-50 


21-26 


energy (keV) 


- 2-4 


- 9-20 


8.0-9.8 


PCA Epoch 4 








channel 


0-5 


16-42 


20-22 


energy (keV) 


- 2-4 


- 9-20 


8.8-9.7 


PCA Epoch 5 








channel 


0-5 


17-44 


21-22 


energy (keV) 


- 2-4 


- 9-20 


9.3 - 9.7 



a Only the upper channel and energy boundaries are given, the 
lower boundary is - 2 keV in all datasets. For epochs 3 to 5, the 
upper boundary varied according with the channel binning used in 
the observation, and the ranges of upper channels and correspond- 
ing energies are given. 



the available data mode). The lightcurves used in determining 
count ratios were extracted with 16 s time resolution. To al- 
low comparison between observations during different epochs 
and using different numbers of PCUs, the lightcurves were all 
normalized to one PCU. 

The data taken from the ASM instrument are the one-day 
average count rates. All three bands have been combined, and 
the energy range is ~ 2 - 12 keV. 

2.2. Calculating the power density spectra 

To calculate the power density spectra the high-resolution light 
curves were divided into segments of 2 13 bins (~ 82 s) and a 
Fourier transform performed on each segment. The resulting 
PDS were then gathered into groups of 10 and each group av- 
eraged into one final PDS. Unfortunately, calculating the PDS 
over short segments results in a low frequency resolution and 
we were therefore not able to model any narrow features in the 
PDS. Reducing the number of averaged segments increases the 
statistical fluctuations. However, our approach allows changes 
in the PDS on timescales as short as ~ 15 minutes to be moni- 
tored. Our preliminary studies showed that such changes do in 
fact occur, and as a result we accepted the lower signal-to-noise 
ratio. This is further discussed in Sect. 13. II 

The frequency range in our observations is 0.01-25 Hz. The 
PDS are corrected for the Poisson level. In principle, dead-time 
effects should be considered when determining the Poisson 
level. For the frequency range used here this correction is very 
small ($ 2%) compared to the statistical fluctuations, and we 
therefore disregard it. As the effect of deadtime increases with 
frequency, we set our upper limit to 25 Hz. For a discussion 
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on dead-time influe nce on PDS see IVikhlinin et all dl994l) and 
IZhang et alJ ill 9951). and for its relevance in the case of the PCA 
detector lZhang et all i ll 9961) and lJernigan et a 

The PDS normalizati on used in this paper is that of 
Belloni & Hasinger ( 1990) and Miyamoto et al. In this 

normalization, integrating over a frequency range gives the 
square of the fractional root mean square (RMS) variability in 
that frequency interval. The PDS are shown both in units of 
power versus frequency and fre quency times power versus fre- 
quency (fPj versus /, see e.g. lBelloni et aljfl997l) . The latter 
method of presenting the PDS more clearly shows the contribu- 
tion of variations per frequency decade to the total RMS of the 
source. The PDS have been rebinned to semi-logarithmically 
spaced bins. 

The observations were provisionally divided into groups, 
roughly depending on the state of the source. Figure ^ shows 
the ASM lightcurve of Cyg X-l, and the vertical lines show 
boundaries between the different groups. The intervals around 
the transitions are deliberately chosen wide in order to encom- 
pass hard and soft state PDS before and after each transition. 
A period of intense flaring activity in the hard state (and pos- 
sibly one or more brief transition s to the sof t state during MJD 
51 830 - MJD 51 960, see ICui et alJEool was treated sepa- 
rately. The reduction resulted in a total of 2016 PDS, of which 
1324 were from the hard state, 226 in flares, 191 during transi- 
tions, and 275 in the soft state. Approximately 2% of the com- 
puted PDS were of too low quality to be used in the analysis. 

2.3. Modeling 

The PDS of Cyg X-l have traditionally been described 
using a double-broken power-law in t he ha rd state (e.g. 
iBelloni & Hasing erl 1 199(1 iNowak et all 1 19991) an d an f' 1 
slope broken at higher frequencies in the s oft state JCui et alJ 
Il997at 119991: IChurazov et all l200ll) . Building on re- 
sults from comparisons between neutr on star (NS) and 
black hole candida te (BHC) sources JPsaltis et alJ |l999; 
IWiinands & van der Klisll 19991 and references therein) it has 
been shown that models using Lorentzian pro files provide a 
good description of the hard state of Cyg X-l dNowakl EoOO. 
hereafter N00: |P03l) . However, these models do not fit the soft 
state of Cyg X-l, and during large flares and state transitions 
an additional power-law component is required 

While there is evidence for several Lorentzian profiles in 
the P DS of Cyg X- 1 , generally only two dominate at a given 
time (N00; Gilfanov et al. 1999). As a base for our model we 
therefore consider the sum of two Lorentzian profiles of the 
form 



Table 2. Combinations of components used. 



W) = 



HiWm 



(/ - v,) 2 + fWm 



(1) 



where Hj is the value of fPf at the peak frequency v,-. In the 
fPf representation, the dimensionless parameter W, gives a 
measure of the width of the profile. The parametrization is 
chosen such that for a given W, and Hi, the integral of the 
Lorentzian (and thus its fractional contribution to the total RMS 
variability) is the same for any value of v,. Note that whil e this 
parametrization differs from others previously used (e.g. N00; 



Model 


Components 


States used 


1 


2 Lorentzians 


hard 


2 


Power-law+2 Lorentzians 


all 


3 


Power-law+Lorentzian 


soft 


4 


Power-law 


soft 



P03), it has the advantage of having a direct interpretation of 
the parameters in the fP j representation. 

We have chosen to use the peak frequency, which is the 
frequency where the contribution of the Lorentzian to the to- 
tal RMS is at its maximum, rather than the centroid frequency, 
since this seems to be the important parameter in terms of fre- 
quency relations van Straaten et al .120021) . In the case of 
narrow features, the peak and centroid frequencies will be vir- 
tually identical. For broad features (such as breaks in the stan- 
dard PDS representation), it is the peak frequency which best 
represents the characteristic frequency of the feature. In such 
cases the centroid frequency has no direct physical meaning, 
and will approach zero. For further discussi on on the peak and 
centroid frequencies see lBelloni et alJ J2002). 

The fitting is done in the fPt parameter space, where our 
parameters have the interpretation described above. Conversion 
to 'standard' Lorentzian parameters shows the advantage of us- 
ing the peak frequency. When the width parameter W, is greater 
than 2, the centroid frequency for the Lorentzian in the PDS 
space becomes negative. Although not a problem mathemati- 
cally, this underscores the fact that in the context of broad com- 
ponents it is the peak frequency that has a physical interpreta- 
tion. In principle we could have frozen the centroid frequency 
to zero in such cases, but we left it free in order to better follow 
the variations of the PDS at very low frequencies. 

In addition to the Lorentzian profiles, when necessary 
(large flares, transitions and in the soft state) we have included 
a cut-off power-law component of the form 



Pl(f)=Af- a e 



(2) 



where A is the normalization constant, a is the power-law 
index, and f c is the turnover frequency. Studies at longer 
times cales of Cyg X-l using RXTE/ASM data dReig et alJ 
2002) show that at low frequencies (below ~ 10 -5 Hz) the PDS 
is approximately described by a ~ / _1 power-law, and there is 
evidence that the PDS follows this index to higher freque ncie s 
during transitions and in the soft state (Zdziar ski et alj E004). 
A simple f~ l power-law does not fit the soft state PDS above 
~ 10 Hz, and we therefore introduce an exponential cut-off. 

As will be described in Sect. [3] when the second Lorentzian 
shifts to higher frequencies in the soft state, its contribution di- 
minishes and it becomes difficult to constrain. Some PDS were 
therefore fit with only one Lorentzian profile plus a power-law 
component. In addition, the canonical soft state PDS exhibiting 
only a / _1 slope do not require any Lorentzian profile. Table|3] 
lists the component combinations used in each state, and exam- 
ple fits for each of the four models are shown in Figs.[2]to|6] 

While it is clear that Cyg X-l exhibits several spectral 
states, there is no unique classification between them. Adding 
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Fig. 1. ASM lightcurve for Cyg X-l from MJD 50 135 (Feb 22, 1996) to MJD 52763 (May 4, 2003). The different groups are 
marked in the plot, with H signifying hard state, Tu and Td transitions from hard to soft and soft to hard respectively, S the soft 
state and F signifies a period of large flares in the hard state. 



to the difficulty are the large flares or 'failed state transitions', 
where the source appears to begin a transition but not complete 
it. We have chosen to call these episodes flares, and reserve 
the use of transition for cases when the source changes state. 
In this study we have attempted to use the PDS when deter- 
mining the state of the source. For the canonical hard and soft 
states, there is no difficulty as their PDS are clearly separate 
and appear only in their respective state. The situation during 
the transitions (and flares) is by its very nature more ambigu- 
ous. In the following, we distinguish between hard and transi- 
tional state PDS by the power-law component appearing only 
in the latter. The same criterion is used to distinguish flares, 
with the difference between flares and transitions being that the 
source does not change state after a flare. The distinction be- 
tween transitional and soft states is not as clear-cut. We have in 
this study assigned the PDS where the second Lorentzian com- 
ponent is not seen or significantly weaker than the first one to 
the soft state, along with the canonical PDS consisting only of 
the power-law component. In our classification, Fig.[2]is a hard 
state PDS, Fig.|3]a PDS from a flare or transition, and Figs.|4]to 



|6|soft state PDS. We shall compare this classification to others 
made from a spectral standpoint in Sect. 13.41 



2.4. Hardness and Spectral Index 

The energy spectrum of Cyg X-l in the hard state can be 
represented as a soft thermal component and a hard power- 
law component with index F ~ 1.6 and a high-energy cut- 
off. There is evidence of an iron fluorescence line at 6.4 
keV and a broad Compton reflection component around 30 
keV (e.g. iBamcinska-Church et aljfl99l lEbisawa et alJll996l 
Gierlihski et al. 119971) . The soft state is dominated by a strong 
thermal component below ~ 3 ke V, combined with a power-law 
continuum extendin g above 1 MeV. The iro n line is visible also 
in the soft state (e.g. Gierlihski et al. Il999i) . The soft radiation 
is expected to come from an optically thick, geometrically thin 
accretion disc, while the hard emission is usually explained by 
a hot, opticall y thin medium where the soft disk photons are 
Comptonized (Shapiro et al. 1976; Poutanen 2001). 
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Fig. 2. Example of fit to hard state PDS using two Lorentzian 
profiles (model 1). This is the most frequent state of the PDS, 
and the majority of our observations are fit with this model. 
Hardness ratios (9 - 20 keV/ 2-4 keV) for this model fall in 
the range of 0.50 — 1.1. The example is from an observation 
made on MJD 50180. 




Fig. 3. Example of fit to a transition PDS using a cut-off power- 
law and two Lorentzian profiles (model 2). The PDS of Cyg X- 
1 take on this appearance during state transitions and larger 
flares in the hard state (failed transitions). Hardness ratios for 
these fittings fall in the range of 0.25 - 0.45. The example is 
from an observation made on MJD 50226. 



We used the hardness ratios (or more specifically count 
ratios) between the 9-20 keV and 2-4 keV channels to pro- 
vide a means of correlating the PDS properties with changes 
in the energy spectrum. To compute a hardness ratio, the 16 s 
lightcurves were first divided according to the start and stop 
times of the data used in each PDS. A mean count rate was 
then calculated in each energy range for each interval, and the 
hardness ratio HR computed. This way, each PDS could be tied 
to an average count rate and hardness ratio. 

Figure shows the flux-hardness relation. In Cyg X-l, 
the correlation between flux and hardness differs between the 
states. During the hard state it is generally negat ive, while in 
the soft state there is a positive correlation dWen et alJl200ll 
IZdziarski et ai1 l2002). There is a distinct peak around HR ~ 
0.3, indicating a rather sharp change in behavior. 

The hardness as defined above is often a useful tool, as it 
provides an easy way to measure spectral changes with very 
good time resolution. However, it cannot readily be compared 
with hardness values from other instruments or energy bands. 
To facilitate comparison with other studies, we derived an em- 
pirical relation between our definition of hardness ratios and 
the spectral index T. This was done by modeling a number 



of energy spectra in the 2-20 keV range with a power-law in 
xspec, as well as calculating their average hardness. While the 
energy spectrum of Cyg X-l, especially in the soft state, cannot 
accurately be described by a simple power-law, this approach 
provides a good approximation of the spectral characteristics of 
the source. The resulting empirical relation between the index 
F and hardness ratio HR was found to be well described by a 
power-law (see Fig.[S}. 

HR = 3.9 xT 2 * (3) 

The transitional state of Cyg X-l has a typical index in 

the range 2.1 ;$ T £ 2.4 i n the energy band 3-12 keV 

dZdziarski et alJl2002l Fig. 3 in lZdziarski & Gierlinskll2004 . 
Comparing to the empirical relation of Eq. [3] this would cor- 
respond to a hardness 0.3 ;S HR S 0.5, with the lower edge 
matching the peak in Fig. The simplicity of the relation in 
Eq.|3]means that spectral index and hardness ratios can be used 
interchangeably, and we will consistently use the latter in this 
paper. 
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Fig. 4. Example of fit to soft state PDS using a cut-off power- 
law and two Lorentzian profiles (model 2). Almost 20 % of the 
soft state observations also allowed a fit with this model. The 
Lorentzian components have shifted compared to Fig. |5J and 
the power-law component extends to higher frequencies than 
during transitions. Hardness ratios for these fittings fall within 
the range of 0.20 - 0.30. The example is from an observation 
made on MJD 52 341. 



3. Results 

In this section we present the results of our analysis. We be- 
gin with a general description of how the PDS evolves during 
spectral transitions, and then describe in detail the relevant pa- 
rameter correlations. 

While narrow, well-defined peaks in the power spectrum 
are generally called QPOs, there is no set definition. This leads 
to a sliding scale where wider features are sometimes referred 
to as breaks and sometimes as QPOs. To avoid confusion and 
connection with the different types of physical processes asso- 
ciated with breaks and QPOs, we will consistently use 'tem- 
poral feature' or 'characteristic frequency' to describe the fea- 
tures seen in the PDS, and in the context of our modeling com- 
ponents refer to them as 'Lorentzians'. For a re view of tem- 
poral feature s in both BHC and NS systems, see Ivan der Klisl 
i 19951 12004 . A possible identification of the different fea- 
tures in BHC in t erms of those found in NS is presented by 
lKlein-WollH2004l> . 



Hz power spectral components in Cygnus X-l 
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Frequency [Hz] 

Fig. 5. Example of fit to soft state PDS using a cut-off power- 
law and one Lorentzian profile (model 3). Almost half the soft 
state observations allowed a fit using this model. Hardness val- 
ues for this model fall in the range of 0. 14 - 0.26. The example 
is from an observation made on MJD 52 341. 

3.1. PDS of the different states 

As Cyg X-l spends most time in the hard state, the majority 
of the observations were made during this state. It is well de- 
scribed by two Lorentzian profiles and no power-law compo- 
nent (see Fig.|2). The two peaks shift in frequency, moving on 
timescales of hours. In very few (^ 1%) of the hard state cases, 
model 1 did not provide a good fit due to a third component 
being evident at higher frequencies. These cases correspond to 
the observations with highest hardness, and will be discussed 
further below. 

During larger flares and transitions, the two peaks shift to 
higher frequencies and it becomes necessary to include the 
power-law component to achieve an acceptable fit. The power- 
law component only significally contributes at lower frequen- 
cies (Fig.|5|l. When analyzing the data we separated flares, tran- 
sitions from hard to soft state, and transitions from soft to hard 
state in an attempt to detect possible differences between them. 
However, in the analysis done in this paper the three appear 
identical. Our resu lts on this point agree with those f ound pre- 
viously re.g.lPottschmidt et alfcOOfltlCiii et all2002tlpol. 

It is during the transitions and flares that the power spec- 
trum of Cyg X-l shows the most rapid changes. We are able 
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Fig. 6. Example of fit to soft state PDS using only a cut-off 
power-law (model 4). This is the canonical soft state PDS 
appearance for Cyg X-l, however only 35% of the observa- 
tions made during the soft state allowed a fit using this model. 
Hardness values for this model fall in the range of 0.1 1 - 0.23. 
The example is from an observation made on MJD 52 257. 



to see changes in the PDS down to our limit of ~ 15 minutes, 
justifying our choice of a short integration time in the analysis. 

In the soft state the power-law component dominates with 
an index close to -1. However, in the majority of the cases 
(~ 65%), one or two Lorentzian components are needed in ad- 
dition to the power-law in order to fit the PDS. The Lorentzian 
components are particularly frequent in the 2001/2002 ex- 
tended soft state. During this period, the source made at leas t 
two brief transitions down to the hard state (also noted by P03). 
No trace of these transitions is seen in the ASM lightcurve of 
Fig. ^ but both the hardness and PDS reverted to that of the 
canonical hard state. 

3.2. Relations between the peak frequencies 

The relation between v\ and V2 is shown in Fig. [9] Our observa- 
tions confirm the correlation between the peak frequencies of 
the PDS found in previous studies of Cyg X-l in the hard state 
(N00; P03). Since our study is based on a much larger data 
set, we can make a more accurate analysis of the hard state 
frequency variations. We are also able to add the behavior of 
the peak frequencies during transitions, flares and into the soft 
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Fig. 7. Flux-hardness correlation of Cygnus X-l. Each point 
represents the mean values of hardness and flux during one 
PDS. Grey points are the hard state, black points the soft state 
and open squares points from transitions and flares. The change 
in correlation around HR ~ 0.3 is clearly visible. The flux is the 
sum of the two channels. 
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Fig. 8. Empirical relation between spectral index F and hard- 
ness HR, obtained by fitting the energy spectra between 2-20 
keV with a power-law. Errors are within the boxes. The equa- 
tion of the power-law is HR — 3.9 x F~ 2 8 . 

states. A clear feature in Fig.|5]is the steepening of the rela- 
tion at V2 - 5 Hz. This corresponds to the region where the 
source goes from hard state into the intermediate/transitional 
state. We note that this is the region where P03 see a shift in 
power from their lower frequency Lorentzian (Li), to a higher 
frequency one (L3). Since we use only two Lorentzians, the 
break could in principle be attributed to a change in the com- 
ponents being fitted. We have however carefully monitored this 
region, and find no evidence of such a change in our data. Since 
our method of analysis allows us to detect changes in the PDS 
on timescales as short as ~ 15 minutes, it is more likely that 
the break is due to a real change in the behavior of the compo- 
nents. The fact that the points lie in a continuous distribution, 
both in Fig.|5]and in the following figures, further supports this 
identification. 

In the soft state, the relation between the two peak frequen- 
cies is somewhat uncertain. We see two major effects behind 
the large spread in these data. Since we do not correct the 
Poisson noise level for dead-time effects, our frequency range 



M. Axelsson et al.: Evolution of the 0.01 - 25 Hz power spectral components in Cygnus X-l 



„ 1 

N 

::: 

£ 0.5 
0.05 




v 2 ( Hz ] 



Fig. 9. Correlation between the peak frequencies of L\ and L2. 
The gray and black points represent the hard and soft states 
respectively, and the open squares are the transition state. The 
break away from the hard state behavior as Cyg X- 1 enters the 
transitional state is clearly visible. 



is limited to the region where these can be neglected. As the 
second Lorentzian approaches our upper frequency limit of 25 
Hz, it becomes difficult to constrain. Even with this rather con- 
servative limit, we cannot rule out that the points near the limit 
are affected by dead-time effects. A second effect which further 
complicates the modeling is the fact that the power of the fea- 
ture is reduced at higher frequencies. As we average over only 
ten power spectra, our signal-to-noise ratio does not allow us to 
model weak components. The spread of points in the soft state 
is therefore due to both the limitations set by dead-time effects 
and the signal-to-noise ratio. 

3.3. Comparing the two Lorentzians 

The strong correlation between the peak frequencies of the two 
Lorentzians in Fig.|9]suggests a close physical link between the 
components. However, when studying the relations between Wt 
and v, the two components show completely different behavior 
(Figs. [TUIandll llrespectivelv). The first Lorentzian shows no 
correlation between these parameters, and W\ remains fairly 
constant. For the second Lorentzian, W2 changes roughly two 
orders of magnitude, and there is a strong correlation, which 
can approximately be described by a power-law with index ~ 
—1.5. In contrast with the relation in Fig.|9j the points for the 
transitional and soft states seem to follow the same relation as 
in the hard state. 

During the hard state there are cases where the Lorentzians 
shift to very low frequencies. These cases correspond to obser- 
vations with very low flux and a high hardness ratio. To provide 
a good fit in these cases, we found it necessary to freeze the ra- 
tio between the peak frequencies by extrapolating the relation 
found at higher frequencies (see Sect. l4.2l below). For these low 
frequencies in the hard state, there is a hint of flattening of the 
index in the W2-V2 correlation (open triangles). Since this oc- 



: 
5 



0.5 



0,1 

0.05 - 




;:-izj 



Fig. 10. Relation between W\ and vi . Symbols are the same 
as in Fig. |9] with open circles being points from the sin- 
gle Lorentzian in model 3. The width is generally constant 
throughout the frequency range. 




Fig. 11. Relation between W2 and v%. Symbols are the same as 
in Fig.|9] with open triangles being points from the hard state 
where the relation between Vi and V2 was fixed. The negative 
correlation between W2 and V2 is clear, and can approximately 
be described by a power-law with index ~ —1.5. 



curs in the region where vi was fixed in relation to V2 we cannot 
rule out that this is an artifact of the fitting procedure. 

In Fig.EB the open circles represent the single Lorentzian 
of model 3 (see Table[2). We identify this component as L\ . The 
identification is supported by the fact that the points from this 
component well matches the other soft state points, both in this 
and subsequent figures. If the component is instead identified as 
L2, the distributions occupy completely separate regions in the 
parameter space. We surmise that as both Lorentzians move up 
in frequency, eventually L2 becomes too weak and moves too 
high in frequency for us to see in our analysis, and only L\ is 
visible. 
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Fig. 12. Power at the peak frequency versus peak frequency for 
L\. Symbols are the same as in Fig.^| The change in behavior 
during and after the transition (open squares) is rather gradual. 



Comparing the relations between Hj and v, shows similar 
behavior for L\ and Li (Figs.ll2landll3>. Both components stay 
at roughly constant height at lower frequencies, and decrease 
rapidly with frequency at higher v,. This change in behavior for 
Li occurs as Cyg X-l enters the transitional state, suggesting 
that the mechanism behind L2 starts to weaken as the transition 
begins, perhaps due to a change in the accretion geometry. We 
surmise that the behavior of Hi at lower freq uencies is most 
likely explained by the effect first noted by Belloni & Hasinger 
(1990) in their analysis of the PDS in terms of a doubly broken 
power-law model. They found a negative correlation between 
the frequency of the first break and the power at that frequency. 
Transforming this correlation into the fPf representation leads 
to approximately constant H\ as a function of vj . Our results 
are also consistent with the relation between normalization and 
pe ak fr equencies of the Lorentzian components seen in Fig. 6 
of Pojl but note the difference in parameters used and identifi- 
cation of components. 

Combining the behavior of //, and W, once again shows 
the difference between the two components. Since both Hi 
and Wi stay roughly constant at low v\ , the contribution from 
Li to the total RMS variability remains the same, weakening 
only above ~ 5 Hz, when the source is in the soft state. In 
contrast, the RMS contribution from L2 falls with V2 over the 
whole frequency range. As H2 starts to decrease, L2 weakens 
very rapidly. Interesting is also that the rapid decrease in height 
starts at //, ~ 5 Hz for both Li and L>2- A possible explanation 
for the similarities and differences in behavior of the compo- 
nents will be presented in Sect. 14. 31 

3.4. Evolution of the PDS 



As described above in Sect. 13.11 there are systematic changes 
in the PDS tied to the changes in hardness. Fi gures IT4l and [TBI 
show the evolution of v, with hardness. Both peak frequencies 
increase with decreasing hardness. This movement follows ap- 
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Fig. 13. Power at the peak frequency versus peak frequency 
for L,2- Symbols are the same as in Fig.^TJ Note the sudden 
change in behavior as the source enters the transitions (open 
squares). This suggests the underlying physical mechanism re- 
acts sharply to changes, e.g. in the accretion geometry, occur- 
ring at these periods. 



proximately the same pattern in all states while the components 
are detectable, which we take as another indication that we are 
following the same two components in all cases. We note that 
this gradual increase in peak frequency would be expected if 
the Lorentzian components are assumed to arise in a region 
that is closer to the central object in the soft state than in the 
hard state (e.g. the inner edge of a truncated accretion disk or a 
transition region betwee n disk and corona). S u ch models have 
been p rese nted bv e.g.lPoutanen etafl dl997t) . IChurazov et alJ 
J200lL and lZdziarski et alJ d2002h . 

Beyond the necessity to include the power-law component 
in our models below a hardness of ~ 0.4, Figs. ^] and ^] 
show no change in behavior around the transition phase hard- 
ness, HR ~ 0.3 seen in Fig. [71 We therefore conclude that the 
changes responsible for altering the flux-hardness correlation 
do not affect the movement of the peak frequencies. The other 
parameters do not show any evolution with hardness, beyond a 
trend for the power-law to extend to higher frequencies as the 
hardness decreases. This is a natural effect of the model as the 
Lorentzian components move away from the lower frequen- 



cies. The behavior is however consistent with the / power- 
law gradually extending hig her in frequency as the source goes 
from hard state to soft state (Zdziarski et al.l2004l) . 

When determining the state of the source, we have done so 
from the point of view of the PDS, as they are the base of our 
analysis. Figures^landElshow mat mere is a unique relation 
between a point in the PDS parameter space and the hardness. 
We are therefore able to compare the 'PDS states' with defini- 
tions based on the spectral index of the source. We see that our 
transitional state PDS roughly fall in the range 0.3 ^ HR ;S 0.5. 
Below this range is the soft state and above the hard state. From 
the relation in Eq. [3] the corresponding range in spectral index 
is 2.1 $ F $ 2.4, with the hard state having a lower and the 
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Fig. 14. Relation between hardness and v\ . The symbols are the 
same as in Fig.^| There is no apparent change in behavior as 
the source enters the transitional or soft states. 
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Fig. 15. Relation between hardness and V2- The symbols are 
the same as in Fig. ^2 Note the tendency for the hardness to 
'flatten out' at ~ 1.1. 



soft state a higher index. We note however, that these bound- 
aries are not strict, and that there is some overlap between the 
diffe rent states. This range matches well with prev i ous re sults, 
e.g. lZdziarski et alJd2002l) : IZdziarski & Gieriiriskjlll2004 . 

To illustrate the complete evolution of the PDS, Fig. [H^ 
presents a series of fittings showing the gradual change from 
normal hard state to canonical soft state PDS. Unfortunately, 
the observations made did not capture the complete evolution 
in a chronological series, so Fig.^]is a montage of PDS from 
observations made during a brief transition to the hard state in 
the extended soft state of 2001/2002. The top three panels are 
from observations made on MJD 52 330 and panels d and e are 
from an observation made on MJD 52 324. The last PDS (panel 
f ) is made on MJD 52 358, when the source had once again re- 
turned to the soft state. Figure^Jillustrates the shift in the peak 



frequencies of the Lorentzians, and the gradual increase of the 
power-law as the source goes from hard state to soft state. 

4. Discussion 

The results of our analysis show that the PDS of Cyg X-l can 
be successfully modeled in all the spectral states using only 
a few components. The observed behavior of the Lorentzian 
components lends weight to the argument that they indeed rep- 
resent the temporal characteristics of the physical processes 
giving rise to the PDS. 

4.1. Identification of components 

As noted in Sect. I3.2I the identification of components in the 
PDS is not trivial. P03 in their study use four Lorentzians to fit 
the power spectrum. While their fourth component, L4, is above 
our upper frequency limit, our model will fit the two dominat- 
ing components at a given time. According to the identification 
of tP03l this means that in the canonical hard state we fit their 
first two Lorentzians, but in the larger flares and transitions 
we fit their second and third Lorentzians. The first Lorentzian 
is then quite weak and will be 'absorbed' in our power-law 
component. Supporting this identification is the change in the 
correlation of the peak frequencies during transitions and large 

flares. 

Recently. | van Straatenl J2004I) has argued that the identi- 
fication of P03| leads to significant deviations during 'failed 
state transitions'. For canonical hard state PDS, the first two 
Lorentzians dominate and follow the relation of Fig. but 
Ivan Straatenl f2004) points out that if the second and third com- 
ponents dominate in transitions and large flares, the peak fre- 
quencies of the first two components deviate significantly from 
the correlation in the normal hard state. If instead the dominat- 
ing components are still identified as the first two Lorentzians, 
the peak frequencies will follow the trend seen in Fig. ^] 
Although the points do not lie along the extension of the power- 
law described at lower frequencies, thei r behavior matches 
that of the other sources, as seen in Wiinands & van der Klij 
( 1999). To support the identification, van Straaten (2004) adds 
data from a flaring episode in November 2000, but a gap be- 
tween the canonical hard state PDS and these flares still re- 
mains. 

Since the study reported here allows for shorter timescales 
to be monitored, we are able to follow changes in the 
Lorentzians in great detail. From Fig. [HO it is clear that our 
model can explain the evolution of the PDS, and provide exten- 
sive coverage of all states of the source. As seen in Fig. II Hour 
results show a continuous distribution between the hard state 
and transitions, and we are able to fill in the gap in the data of 
Ivan Straatenl d2004l) . The r esults here therefore favor this iden- 
tification over that of |P03l 

4.2. Comparisons with models 

The temporal features observed here have been detected both 
in NS and BHC systems. The origin of these characteristic 
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frequencies has been long debated in the literature, and vari- 
ous models have been pr oposed in eac h physical scenario (for 



0.01 0.05 0.1 0.5 1 5 10 

Frequency [Hz] 

Fig. 16. Montage showing the evolution of the PDS of Cyg X-l 
from hard state (panel a) to canonical soft state (panel f). The 
top three panels (a-c) are from observations made within a few 
hours during MJD 52 330, and clearly illustrate the changes 
that can occur on this timescale. Panels d and e are from MJD 
52 324, and are also no more than a few hours apart. The bot- 
tom panel (f) is from an observation on MJD 52 358. Note the 
weakening and narrowing of the upper frequency Lorentzian 
(dashed line) as it shifts to higher frequencies. 



an extensive review, see Ivan der Klis 2004). However, recent 
comparative studies between these two types of sources have 
shown that they present very similar freq uency correlations, 
over a frequency r ange of many decades (Psa ltis et alj Tl999: 
Bello ni et al . 2002). The similarities suggest the same underly- 
ing physical mechanism, which cannot be dependent on either a 
magnetosphere, a solid surface, or an event horizon dWiinandsl 
Beat-frequency models, which have been extensively 
used in the context of NS, can not easily be extended to black 
holes as these are unable to provide an anchor to the magnetic 
field, a condition that is needed to generate a beating with the 
Keplerian frequency. 

Several classes of models have been pro posed to e xplain the 
observed QPO's and their correlations (see lPsaltis et al.lfl999l 
and references therein). They commonly ascribe the variations 
to processes in the accretion disk. This can be done by identi- 
fying observed frequencies with orbital, epicyclic and preces- 
sional frequencies at certain radii (such as the inner edge of the 
disk). Another approach is to study disk oscillation modes. For 
example, in the region close to the compact object, global disk 
oscillation modes can become trapped in particular di sk annuli. 
For co mpre hensive summ aries on this topic, se e e.g.lWagonerj 
( | l999j) and iKatol d200ll) . As demonstrated in Wa goner et alJ 
(2001J), disk mode frequencies are mainly determined by the 
mass and angular momentum of the compact object, and de- 
pend only weakly on the disk parameters. While this property 
makes the models good candidates in ex plaining the high co- 
herence seen in higher frequency QPO's ( Barre tet al.l2004l) . it 
also makes it difficult to explain the large (over one order of 
magnitude) shifts in peak frequency observed in our study. In 
addition, the features we model with Lorentzians have low co- 
herence. Our parametrization does not directly give a value for 
the more generally quoted quality factor Q = v/Av, but for L\ 
this value is S 1 in all states. For Li, Q varies with the spectral 
state. It is low (Q < 1) in the hard state and while higher in the 
soft state, Q < 5 in all observations. 

In the following, we will consider the somewha t idealized 
case of the r elativistic precession model (RPM, Stell a & Vietril 
M 9981 M 9991) . It states that the timing features are the observ- 
able effect of the relativistic nodal and periastron precessions 
of a disk element or a 'blob' in a slightly tilted and eccen- 
tric orbit, respectively. Equally possible is the interpretation 
of the features arising in a narrow band o r annulus in the in - 
ner region of the accretion disk (see e.g. Ivan der Klisl Eo04). 
Although the predicted frequencies are those for free particle 
orbits, thes e can still often be reco gnized in the models of disk 
modes (e.g. Wagoner 1999). Psaltis & Norman (2000) showed 
that a sharp transition region in the accretion disk acts as a 
low-band pass filter with strong resonances near the frequen- 
cies predicted by the RPM, motivating a closer look at this 
model. IStella. Vietri. & Morsinkl ill 9991 hereafter SVM) have 
shown good agreement between the predicted frequency corre- 
lations in the RPM model and a small number of BHC and 
NS observations. Although the model currently requires NS 
spin rates and masses higher than measured, these problems 
are dependent on the chosen equation of state, and the RPM 
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model successfully predicts the observed quadratic dependen- 
cies betwee n lower and higher frequency features seen in many 
sources (see lvan der Klisl2004 and references therein). Based 
on previous calcu lations (bv lBardeen et alJl972llOkazaki et alJ 
ll987tlKatol 19901) of the Keplerian epicy clic v r , and vertical 
vg frequencies for a Kerr black hole, SVM derived that in the 
weak field approximation for a slowly rotating black hole 

Vnod oc fl,M" 5 V^ , (4) 

where v no d = - vg is the nodal precession frequency, v per = 
- v r is the periastron precession frequency, M is the mass, 
andfl, = a/R g is the dimensionless specific angular momentum 
(with a = J/Mc, R g = GM/c 2 being the gravitational radius, 
J the angular momentum, c the speed of light in vacuum, and 
G the universal gravitational constant). We used a power-law 
model to fit the correlation between frequencies for the hard 
state data (see Fig. 1171. up to the break observed in the transi- 
tional state at V2 - 5 Hz. The best-fit parameters did not change 
significantly selecting a somewhat smaller, more conservative 
upper limit, nor was the fit affected by ignoring the points at the 
lowest frequencies. We obtained a power-law index 1 .20 + 0.01 
in remarkably close agreement with Eq.0] 

In lSVMl it is argued that due to the geometry of a tilted inner 
disk, the second harmonic of the nodal precession frequency 
gives the strongest modulation, and that this assumption pro- 
vides the best match for their observations. We also found that 
a good fit of our data using the v\ = v no d identification would 
only be possible for an unrealistically large black hole mass. 
Adopting the IS VMl identification, i.e. in our case v\ = 2v no d 
and V2 = Vper, and using the full relativistic expressions for the 
precessional frequencies, it is possible to derive some of the 
physical parameters of the syste m. For this purpose we used 
the M = 8Mq mass estimate (see lNowak et al J fl 999 . and ref- 
erences therein). Although this M value is somewhat uncertain, 
the frequency relation depends weakly on the black hole mass 
as shown by Eq.|4] At this frequency range, there are two possi- 
ble solutions for the specific angular momentum a*, depending 
on whether the orbital motion is prograde or retrograde. We ob- 
tained a, = +0.49 + 0.01 and a, = -0.57 + 0.01 respectively, 
comparable to previous results based on modeling of energ y 
spectra (see e.g. lShapiro & Lightmanll976tlZhang et alJl997l) . 
In Fig.^Jthe best fit for the prograde case is shown (solid line), 
and some derived orbital radii are indicated. During the hard 
state, within the assumptions of the RPM, the inner radius R m 
of the disk ranges from 50 to 22 R g when reaching the transi- 
tional zone. If the retrograde case is instead assumed the corre- 
sponding values change very little, with the inner disk ranging 
from 55 to 25 R g . These ranges are in very g ood agreement with 
coro nal model estimates for Cyg X-l (see di Matteo & Psaltis 
1999, and references therein). 

The behavior during the transitional and soft states seems 
to be more complex, in part because of the spread of the data 
points at the highest frequencies. If all the data above the 
break is modeled with a second power-law (thus given an index 
1.7 + 0.03) the fit residuals clearly do not appear randomly dis- 
tributed. Another possible interpretation of the data is that after 
the transition, the points corresponding to the soft state follow a 
power-law with an index identical to that of the hard state. This 



picture is reinforced when comparing our Cyg X-l RXTE data 
to those presente d in a study of the frequency co rrelations in 
X-ray binaries by Wiin ands & van der Klis ( 1999), from many 
different NS and BHC sources (open squares in Fig. ll7t . The 
data points seem to lie along 'parallel tracks' as predicted by 
the RPM. This behavior, along with the gradual shift between 
the two tracks during the transitional state, must then be ex- 
plained. It seems likely that a reconfiguration of the accretion 
disk geometry during the state transitions is responsible for the 
changes in the relation between the frequencies. 

A suggestio n to explain the state transi tions in Cyg X- 
1, first made by [Shapiro & Lightmanl (1976) and later on by 
Zhan g et alJ Il997t) . is that they are produced by the reversal 
of the disk rotation. This would occur due to the unstable na- 
ture of the stellar wind from the primary star feeding the ac- 
cretion. In this scenario, the system would switch back and 
forth from a retrograde rotation in the hard/low state, to a pro- 
grade rotation in the soft/high state. Differences in luminos- 
ity would be explained by the changes in the 'last stable orbit' 
/?i ast , that in the case of maximal rotation would change from 
^iast(a* = -1) = 9R g to ^ia S t(fl* = +1) = R g . According to our 
previous fit of the hard state data, a reversal of the disk rotation 
would imply that a* = +0.57 in th e soft state. However, this 
value does not fit the soft state data. SVM showed that correc- 
tions for non-negligible eccentricities or tilt angles are gener- 
ally small. Very large eccentricities (e > 0.85) would need to 
be invoked to explain the data. The rotation reversal scenario 
could thus only be maintained with additional assumptions. 

Reducing the spread of points in the soft state data could 
help determine the behavior of the relation and thereby further 
constrain models. A more thorough study of the soft state data, 
extending to higher frequencies, is currently underway. This 
analysis will also make use of data recently made public from 
an additional soft state in June/July 2003. 

4.3. Behavior of components 

As noted bv lBelloni etaD i2002b . when using Lorentzian com- 
ponents to model both broad and narrow features in the power 
spectrum, the peak frequency v, does not have an immediate 
interpretation in the picture of a damped harmonic oscillator. 
Interpreting the features as arising by superpositions of a range 
of frequencies (such as from an annulus in the innermost re- 
gion of the accretion disk), v, instead becomes a measure for 
the highest frequency contributing, or the one contributing the 
most power. Our width parameter W, can then be seen more 
as a measure of the extent of the contributing region, than as a 
measure of the lifetime of an ind ividual process . For a deeper 
discussion on this topic, see iBelloni et alJd2002l> . 

The different behaviors of the width parameter W, in the 
two components directly rule out the simple case of the compo- 
nents originating in the same 'blob'. If the Lorentzians instead 
arise in an annulus close to the inner radius of the accretion 
disk, increasing v, indicate decreasing R[ n . The decrease of W2 
with frequency may then be understood by a more narrow an- 
nulus dominating the radiation at smaller R ln . However, W\ is 
roughly constant as Vi increases. If the two components arise 
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Fig. 17. Relatio n between peak frequencies (black points) together with the fit to the rel ation predicted in the relativistic p reces- 
sion model of IStella & Vietril dl998l Il999l) . Also plotted are the points from Fig. 2a in lWiinands & van der Klisl (fl999) (open 
boxes), with data from both neutron star (atoll sources) and BHC sources. The arrows mark different inner radii of the accre- 
tion disk, in units of gravitational radii. The solid line is the fit to the 1230 hard state points in the case of prograde rotation, 
when v\ = 2v noc |. The best-fit power-law has an index of 1.20 + 0.01. As the source enters the transitional state, the relation 
clearly breaks away from that predicted by the RPM model. The dashed line corresponds to a power-law of index 1 .20 shifted to 
approximately follow the points from the soft state. 



in the same annulus, W\ may be expected to behave in the same 

way as Wo- As can be seen in Figs.^Sand^J ^ s ' s not ^ e 
case. Any model in which the two components arise from the 
same annulus then requires a mechanism allowing L\ to show 
different behavior than Li. In the context of the RPM, where 
L\ is identified as the nodal frequency, a possible explanation 
is that warping of the disk leads to a smaller deviation from 
the orbital plane at smaller R m . This would lead to the modula- 
tion decreasing with decreasing R\ n , thus compensating for the 
dominance of the innermost ring in terms of emitted radiation. 
The result would be a larger 'effective' annulus contributing to 
L\ than to Li, which is dominated by the narrow innermost re- 
gion. For L\, the net effect must be an approximately constant 
width as seen in Fig. ^2 We note that during the hard state, 
R m decreases by about a factor of 2 (Fig. 1171 . while the flux 
increases with approximately the same factor (Fig.|7J, suggest- 
ing a ~ 1 //?;„ relation. If the deformation of the disk scales 
as R, the net effect would be consistent with a constant power 
of the component as seen for L\. If this simple picture is true, 
L\ could provide important constraints on the geometry of the 
accretion disk close to R m . 

The argument used for the different behaviors of W\ and 
W2 can also be extended to the case of Hi during transitions 
(Figs.[ElandE}- ^2 starts decreasing rapidly as the transition 



period starts, whereas Hi shows a more gradual reaction. If the 
transitional stage is considered to be a phase of rapid changes 
in the accretion geometry, a feature arising from a narrow 'ring' 
will react more sharply than one originating in a wider disk an- 
nulus. After the transition both Lorentzians show the same be- 
havior, indicating that the geometry once again becomes more 
stable. Since the normalization used gives the relative power of 
different components, it is natural that as the Lorentzians shift 
to higher frequencies, the power-law component will start to 
dominate the PDS. The decrease in power of the Lorentzians 
will then mainly be a decrease in relative strength, rather than 
an intrinsic weakening of the physical processes giving rise to 
the components. We note that this scenario is also a possible ex- 
planation for the changes in the frequency correlation in Fig. 1171 
during transitions. During the transition the changes in accre- 
tion geometry disturb the relation between the peak frequen- 
cies, but as the source settles into the soft state the correlation 
is resumed. An explanation must still be found for the differ- 
ence between hard and soft state, e.g. large deviations from a 
circular orbit or a shift from the second to higher harmonics. 

The decrease of W2 with peak frequency is also interesting 
when comparing na rrow QPOs with br oader features (such as 
breaks or 'humps' ). |Psaltis et alJ i 19991) observe that whereas 
so-called lower kHz QPOs seen in neutron star sources are 
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narrow well-defined peaks, components identified as the lower 
kHz QPO at much lower frequencies are broader. As our obser- 
vations of Li indicate, studies such as this one have the poten- 
tial to show the complete evolution from the broader features at 
low frequencies to narrow QPOs at high frequencies. Tracking 
these components will determine whether they indeed originate 
from the same process. 

A feature seen especially in Fig.^]is the tendency for the 
hardness to 'flatten out' at a value HR ~ 1.1. Cyg X-l only 
reaches these levels of hardness in observations with very low 
flux. The behavior of the hardness at these times could be ex- 
plained by the index of the power-law component in the energy 
spectrum (T) reaching a minimum value, and changes in flux 
(and thereby the PDS) arising from another component, e.g. 
small flares. If the variations of F are attributed to variations 
in the accretion disk rather than the hot inner flow/corona (e.g. 
IZdziarski et al.l20o3) . the 'flattening' of the hardness could be 
seen when these variation reached some outer boundary, and 
changes in flux were instead due to smaller variations in the 
processing of the soft input photons. That the peak frequencies 
of the Lorentzian components also vary at these times might 
suggest they are related to the corona rather than the accre- 
tion disk. Arguments for relating 'breaks' or QPO features in 
Cyg X-l to a transition region at the inner e dge of the disk 
or to t he c orona have been presented by e.g. Chur azov et alJ 
il200ll) and P03. However, we note that this behavior can also 
be explained by variations in the loca l disk accretion rate. Such 
a mechanism has been suggested bv IZdziarski et al J d2002l) to 
explain variations in amplitude of the hard state energy spec- 
trum with no change in shape, i.e. roughly constant T. 

As stated in Sect. 12.31 during the observations with lowest 
flux we froze the relation between V\ and V2 according to the 
solid line in Fig.^] The flux of these observations were typi- 
cally below ~ 500 counts/s in the 2-9 keV range, and the hard- 
ness was HR ^ 1 . The vast majority of these observations were 
made between MJD 50400 and MJD 51 000 (November 1996 
- July 1998, cf. Fig.[TJ, and constitute ^ 10% of the hard state 
observations. During these low flux observations, we could also 
see an excess above 10 Hz in almost half the cases. This may 
well be a third component, which we only see the times it has 
moved to the lower edge of its frequency range. In observations 
with a more typical flux level it moves too high in frequency to 
be detected in our analysis. We su rmise this is the highest fre- 
quency component seen by N00, and L\ in |P03l This feature 
could add to the spread of points at low v,, but does not other- 
wise affect our results. 



5. Conclusions 

After more than eight years of observations, the RXTE data 
archive is an excellent source for long-term studies. Using data 
from this archive we were able to follow the evolution of the 
PDS of Cyg X-l from February 1996 to May 2003. By choos- 
ing a simple yet flexible approach for our analysis we could 
model the major features in all states with only three com- 
ponents. Our method also allowed us to detect changes in the 
power spectrum on timescales as short as ~ 15 minutes. 



Our analysis provides several new results. The main result 
is the successful modeling of the PDS of Cyg X-l through all 
the states. We prove that our choice of models and parameters 
successfully enables us to show a natural transition from hard 
state to intermediate to soft state and back, and also allows us 
to study the behavior of the Lorentzian components in detail. 
The combination highlights new correlations between the pa- 
rameters with interesting physical interpretations. 

We find that the spectral state of the source can be uniquely 
determined from the parameters of the Lorentzian components. 
The parameters used here are the width W,- and peak frequency 
V, of the Lorentzians, along with the power at the peak fre- 
quency, Hi. The correlations between parameters can all be de- 
scribed by continuous functions, which differ between compo- 
nents. This is most evident when comparing Figs. II 1 1 and 1 1 01 
While the correlations between the peak frequencies and be- 
tween Hj and v, change with the spectral state (cf. Figs.l9"llT4l 
and!15>. the other correlations show no such behavior. 

A second result of our study was to show that the PDS 
of Cyg X-l is dominated by the same two Lorentzian com- 
ponents at all times. We base this conclusion on the fact that 
we saw no change in the identification of components on our 
analysis timescale and that the points in Figs.|9]to[2]each de- 
scribe a continuous distribution. The relation of the peak fre- 
quencies of these components follows the pattern seen in both 
other BHC and NS systems, and we therefore conclude that the 
physical processes responsible for them cannot be explained by 
invoking magnetic fields, a solid surf ace, or an event horizon . 
A possible model is the RPM of lStella & Vietrll dl998l fl999b . 
which describes the hard state data remarkably well, but re- 
quires some rather ad hoc assumptions to fit the soft state. One 
possibility is to invoke a change of rotational direction during 
the state transitions. We find that this scenario alone is insuf- 
ficient to explain the data, and additional assumptions are still 
required. Another point that speaks against the rotational rever- 
sal model for state transitions is the observed similarity with 
other sources, some of which are believed to accrete through 
Roche lobe overflow and thereby not experience any changes 
in the rotation of the disk. 

We also compared the results of transitions from hard to 
soft, soft to hard and flares for systematic differences. Within 
our analysis we were unable to differentiate between the three, 
and we cannot show any way to predict whether such an event 
will lead to a full source transition or a large flare. 

With this study we have shown the importance of look- 
ing at the PDS of Cyg X-l in a broader context covering all 
the states of the source. To this end we are currently studying 
the soft state PDS in more detail and adding newly available 
data. Determining the behavior of the Lorentzian components 
at higher frequencies will be of great importance for increasing 
our understanding of the physical processes behind them. 
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